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Abstract 

We examine the dependence structure of finite block-maxima of multivariate 
distributions. We provide a closed form expression for the copula density of the 
vector of the block-maxima. Further, we show how partial derivatives of three- 
dimensional vine copulas can be obtained by only one-dimensional integration. 
Combining these results allows the numerical treatment of the block-maxima of 
any three-dimensional vine copula for finite block-sizes. We look at certain vine 
copula specifications and examine how the density of the block-maxima behaves 
for different block-sizes. Additionally, a real data example from hydrology is 
considered. In extreme-value theory for multivariate normal distributions, a cer¬ 
tain scaling of each variable and the correlation matrix is necessary to obtain a 
non-trivial limiting distribution when the block-size goes to infinity. This scaling 
is applied to different three-dimensional vine copula specifications. 

Keywords: Multivariate copula, vine copulas, finite block-maxima, scaled block- 
maxima, extreme-value scaling. 


1 Copula Density of the Distribution of Block-Maxima 


Basically, block-maxima have been used in extreme-value theory as one app roach 


to de rive the family of General Extreme-Value (GEV) distributions flMcNeil et al. 


( 201C)h ). In the recent past the block-maxima method has been studied more thor oughly 
and compared to the peaks-ove r- threshold (POT ) method in iFerreira et al.l (120141) 


and l.laruskova and Hanekl (120061) . iDombrvI (120131) justifies the usage of a maximum- 


likelihood estimator for the extreme-value index within the block-maxima framework. 
The numerical convergence of the block-maxima approach to the GEV distribution is 
examined in Faranda et al.l ( 2011 b Moreover, the block - maxim a method has found 
its way into many application areas: Martv and Blanchetl ( 2012 1 investigate long-term 
changes in annual maximum snow depth and snowfall in Switzerland. Tem pera ture, 
precipitation, wind extremes over Europe are an a lyzed i n iNikulin et al.l (120111 ). A 
spatial application can be found in Naveau et al. (2009). Roccol ( 2014 1 provides an 
overview over the concepts of extreme-value the ory being used in finance. While many 
of the articles use univariate concepts, B ucher and Segers (120131 1 treats how to esti¬ 
mate extreme-value copulas based on block-maxima of a multivariate stationary time 
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series. In contrast to the existing literature (known to the authors), in the following 
we will consider finite block-maxima of multivariate random variables focusing on the 
dependence structure. 


Let U = (U \,..., U d )' be a random vector with 1A\ 0, l]-distributed margins, copula 
C and copula density c. We consider n i.i.d. copies U* = (U^ i,..., U^d)' of U, i = 
1,... ,7i. We apply the inverse probability integral transform to each component of U* 
to obtain marginally normalized data (called z-scale): 

Z i = (Z itl ,..., Z i>d )' with Z id := $ _1 ( U\j) ~ W(0,1), 

for i — 1,..., n, j — 1,..., d, where is the quantile function of the standard normal 
distribution W(0,1). We consider this normalized scale since later we want to compare 
this to the limiting approach used to derive the multivariate Hiisler-Reiss extreme- 
value copula. We are interested in the distribution F^ of the vector of componentwise 
block-maxima 


M (n) = (M[ n \ ..., M< n) )' with ilf (n) := max Z id 


for j — 1,..., d. According to Sklar (119591 ) the dependency structure is determined by 
the corresponding copula C M (»). Since Zij, i = 1 are i.i.d., we know that the 

marginal distribution functions of M- n ' ) are given by 


p ( z i J Z nJ < rrij) = 


and hence the corresponding densities are 

c(n). 


fj ( m j ) = n$(mj) n i p(rrij) 


(1.1) 

( 1 . 2 ) 


for rrij G M, j = 1 Here $ and ip denote the distribution function and the 

density of the standard normal distribution, respectively. Thus, the copula C M (n) is 
the distribution function of 

V = (Hi,..., V d y with Vj := $ (iwf ~ U[ 0,1]. 

For n G N the copula of the componentwise maxima C M („) can be expressed in terms 
of the underlying copula C as follows 

C M (n)K, ■ ■ ■ ,u d ) = c (u{ /n , ... ,«y n ) , (1.3) 

where Uj G [0,1], j = 1 Since C is assumed to have a density c, Equation 

11.31 yields that C M (,o also has a density, denoted by c M (n). Using Sklar’s Theorem, 
Equations II .11 and 11.31 imply that the joint distribution function of is given by 

F (n) (?ui,.. .,m d ) = U M (n) ^(mi),.. .,F] n) (m d )J = ... ,$(m d )) n . 

Theorem 1.1. The density of the copula of the vector of block-maxima satisfies for 
Uj G [0,1], j = l...,d: 


%(«) (wi, • • •, Ud) 
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where d A n := min {d, n}, Sd,j ■— {V\V partition of {1,..., d} with \V\ — j} and 


d M C 




d p C(v i, ...,v d ) 
dv mi • • • dv mp 


1 In 1/n 

Vl =u{ ,...,V d =U d ' 


for M = {mi,.. .,m p } C {1,.. .,d}. 

The proof of Theorem 11.11 as well as all other proofs can be found in the appendix 
at the end of this chapter (pages [14] ff.). 

For the joint density f n > of the block-maxima with marginally normalized data (on 
the z-scale) we also obtain an explicit expression. 

Corollary 1.2. For mj G R, j — 1,..., d, we have 


dAn 


f (n \mi, ...,m d )= ( JJ 'p(rrij) ) ■ ^ 

VJ=1 


n\ 


— •C'(<f>(mi),...,<F(m d )) 


n-j 


{(n-j)\ 

e n d M C($(mi),...,$(m d ))\. (1.5) 

VeS dj Mer j 


Example 1.3. Let d = 3, n G N, i.e. Uj = {U i>u U h2 , U i>3 )', Z, = (Z iA , Z i>2 , Z ij3 )' and 
M( n ) = (M[ n \ M^\ M^ n) y, i — 1,..., n. If n > 3 the copula density of the vector of 
the block-maxima is given by 


c M (n){ui,u 2 : u 3 ) = 


(u\U 2 uf) 


i-1 




I 1/n 1/n 1/n 

nC u{ , u 2 , u 3 ' 


n— 1 


. 1/n 1/n 1/n 

C[U{ ,u{ ,u 3 ' 


n —2 


+ n(n — 1 )C [u / n , u 2 n ■, u l 3 n ^ 

■ \c (u\ /n , uT, Pp) d 23 C yp) 


+ d 2 c („p\ u'p, uT) d a c («;/”, uT, d'" 

+ d 3 C (u\' n , uT, uA d l2 C ( u\ /n , uT 


+ n(n — l)(n — 2 )C (up, up, up) d\C (up, 
■ d 2 c («p,up, up) a 3 c (up, up, up) 


1/n 1/n 1/n 


u 2 , u 3 


2 Vine Copulas 


While the catalog of bivariate copula families (see for example Joe ( 19971) ) is large 
this is not the case for multivariate copula families. They were initially dominated 
by Archimedean and elliptical copulas, however for complex dependency patterns such 
as asym metric dependence in the t ai ls these classes are insufficien t. The class of vine 


copulas ([Bedford and Cookel (1200211 . iKurowicka and Cookel (120061) . lAas et al.l (1200911. 
Kurowicka and Joe ( 201 ll) J can accom modate such patte rns. See Stober and Czado 
( 20121 ) for a tutorial introduction and Czado (,2010) and Czado et al. (12013 ) for re¬ 
cent reviews. Basically vine copulas are constructed using bivariate copulas called pair 
copulas as building blocks which are combined to form multivariate copulas using con¬ 
ditioning. The pair copulas represent the copula associated with bivariate conditional 
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distributions. The conditioning variables are determined with the help of a sequence of 
linked trees called the vine structure. Further it is commonly assumed th at the condi¬ 
tioning value does not influence the copula and its parameter. See Stober et ah ( 2013f ) 


for a discussion of this simplifying condition. Further they show that multivariate Clay¬ 
ton copula is the only Archimedean copula which can be represented as a vine copula, 
while the multivariate t-copula is the only scale elliptical one. For the multivariate 
Gaussian copula and the t-copula the needed pair copulas are bivariate Gaussian or t- 
copulas, respectively. The corresponding parameters are given by (partial) correlation 
parameters. 

Vine copulas allow for product expressions of the density. We only consider three- 
dimensional vine copulas which can be expressed as 


c(«i, U2, U 3 ) — Cu(ui, U2)C2-i{u2 : M 3 )ci 3;2 (Ci | 2 (u\ \ U2) , C 3 | 2 (u 3 1 ^2) ) • 


( 2 , 6 ) 


%3 -k denotes the bivariate copula density corresponding to bivariate distribution 


Here c. 

( Ui , Uj) given U k = u k and 


In,) denotes the conditional distribution function of 


Ui given Uj = Uj, which can be expressed as 


d 

Ci\j(ui\uj) Gy (iii , Uj ) 


Further we write the bivariate copula densities in terms of their copula, i.e. Cij(ui, Uj) = 
dSdu %)■ For the pair copulas Ci 2 , (7 23 , Ci 3;2 arbitrary bivariate copulas can be 

utilized. Many bi variate families including rotations are implemented in the R library 
VineCopula (see Schepsmeier et ah ( 2014h ). which allows for parameter estimation and 
model selection of vine copulas in arbitrary dimensions. 

Now we will consider the three-dimensional case and derive expressions for the 
partial derivatives needed in Theorem 11.11 for the expression of the copula density for 
the block-maxima. 


Theorem 2.1. For the vine copula density (|2.6[) we have: 

1. C(u 1 ,u 2 ,u 3 ) = f 0 U2 C 13;2 (Gi| 2 (ui|u 2 ),G 3 | 2 (ii 3 |u 2 )) dv 2 , 

2. (a) d 1 C(ui,u 2 ,u 3 ) = d 1 C 13 \ 2 (C 1 \ 2 (ui\v 2 ),C 3 \ 2 (u 3 \v 2 ))ci 2 (u 1 ,V 2 )dv 2 , 

(b) 5 2 C(rii, u 2 , %) = Gi 3;2 (Gi| 2 (mi|u 2 ), C' 3 | 2 (u 3 |u 2 )), 

(c) d 3 C(ui,u 2 ,u 3 ) = f Q U2 9 3 Ci 3 | 2 (Gi| 2 (ui|u 2 ), G 3 | 2 (m 3 |u 2 ))c 23 (u 2 , u 3 )dv 2 , 

3. (a) di 2 C(ui,U 2 ,u 3 ) = 9iGi 3;2 (C'i| 2 ('Ui|u 2 ), G 3 | 2 (u 3 |-u 2 ))ci 2 (rii, m 2 ), 

(b) d 13 C(u 1 ,u 2 ,u 3 ) = / Q U2 ci 3;2 (Ci| 2 (mi|u 2 ), C' 3 | 2 (rt 3 |u 2 ))c 23 (u 2 , u 3 )ci 2 (ui, v 2 )dv 2 , 

(c) d2 3 C{Ui 1 U2,U 3 ) = 9 3 Gi 3;2 (Ci| 2 (Mi|u 2 ), G 3 | 2 (u 3 |m 2 ))c 23 (-U 2 , u 3 ), 

4- c(ui,U2,u 3 ) = Ci 2 (tii, ■u 2 )c 23 ('U 2 , n 3 )ci 3;2 (Ci| 2 ('Ui|ri 2 ), C' 3 | 2 ('U 3 |m 2 )). 

Theorem 12.11 shows that the copula density corresponding to the three-dimensional 
vector of block-maxima based on an arbitrary vine copula is numerically tractable 
since only one-dimensional integration is needed. In particular this allows a numerical 
treatment for the block-size n in a finite setting. Additionally we can use the vine 
decomposition for a three-dimensional Gaussian or t-copula instead of requiring three- 
dimensional integration to calculate the corresponding density of the block-maxima. 
Two examples illustrate this way of proceeding. 
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Example 2.2. As a first example we take a three-dimensional Clayton-vine, i.e. all 
three pair-copulas are bivariate Clayton copulas. As parameters we choose S 12 = 6 , 
c ) 2 3 = 7.09 and f ) 13 . 2 = 4.67 corresponding to Kendall’s r values of r i2 = 0.75 and 
r 2 3 = 0.78, 7 "i 3 ;2 = 0.70. Figure 0 shows the copula density of the block-maxima of this 


z 3 = < t > '(0.2) z 3 = <D 1 (0.5) z 3 = <I> 1 (0.8) 



-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 




Figure 1: Two-dimensional slices of the copula density for block-maxima of a three- 
dimensional Clayton-vine with normalized margins (ri 2 = 0.75, r 23 = 0.78, ti 3;2 = 0.7). 

vine with normalized margins (i.e. on the z-level) for block-sizes n = 10,50,10 3 . Each 
row represents one block-size and contains three contour plots. Since it is difficult to plot 
three-dimensional objects in a simple way we decided not to show the isosurfaces but cut 
the three-dimensional object into three slices parallel to the z^-z^-plane. Each column 
presents the contourplot of one slice, where the z^-value is fixed to $ _ 1 (0.2) ; $ _ 1 (0.5) 
or <3> - 1 (0.8), respectively. Furthermore, we plotted the contours of the independence 
copula with normalized margins. One can see that already for n = 10 3 the contours 
of the copula density of the block-maxima with normalized margins practically coincide 
with the ones of the independence copula. 

Remark 2.3. Even though it is not known whether all Clayton-vines lie in the domain 
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of attraction of the independence copula, one can show that the Clayton-copul a, which 


can be represented as a Clayton-vine with specific parameter restrictions tStober et al. 
1\2013()). lies in the do main of attraction of the independence copula. According to 


Gudendorf and Seger\ \201(\) an Archimedean copula with generator p lies in the do¬ 
main of attraction of the Gumbel-copula with parameter 6 := — lim^o E [1, oo) 

if the limit exists. For the Clayton-copula this limit is equal to 1. Therefore, the copula 
of the block-maxima of a Clayton-copula converges to the Gumbel-copula with 6 = 1, 
which is the independence copula. 


z 3 = < t > '(0.2) z 3 = <t>^(0.5) z 3 = <t>~ 1 (0.8) 





Figure 2: Two-dimensional slices of the copula density for block-maxima of a three- 
dimensional Gaussian vine with normalized margins (ri 2 = 0.5, T 23 = 0.57, ri 3 ; 2 = 0.35). 


Example 2.4. The second example we present is a three-dimensional Gaussian vine, 
i.e. all three pair-copulas are bivariate Gaussian copulas. As parameters we choose 
P12 = 0.71, P23 = 0.78 and pi 3;2 = 0.52 corresponding to Kendall’s r values ofr \2 = 0.50 
and t 2 3 = 0.57, Ti 3 ; 2 = 0.35. 
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Figure [| shows the copula density of the block-maxima of this vine with normal¬ 
ized margins (i.e. on the z-level) for block-sizes n = 10 , 50,10 3 . As above each row 
represents one block-size and contains three contour plots corresponding to z 3 -values 
fixed to <f> _ 1 ( 0 . 2 ), 1> _ 1 (0.5) or 1> _ 1 (0.8), respectively. Again we detec t c onverg e nce to 
the independence copula. This is also what one would expect: Hii sler and Reisi (1989 ) 
showed that in order to achieve that the distribution of the maxima of a multivariate 
Gaussian distribution converges to a non-trivial limiting distribution, a proper scaling 
of the margins and the correlation coefficients is necessary. This will be discussed in 
Section 0 


Example 2.5. Hydrology is one of the areas where block-maxima are important. Espe¬ 
cially, the water levels of rivers can be interesting when it comes to analyzing the risk of 
floods. We consider a three-dimensional data set containing the water levels of rivers 
in and around Munich, Germany, from August 1, 2007 to July 31, 2013. The data 
has been taken from Bavarian Hydrological Service (http://www.gkd.bayern.de). 
The three variables denote the differences of the 12 hour average water levels at the 
following three measuring points: the Isar measured in Munich, the Isar measured in 
Baierbrunn (south of Munich) and the Schwabinger Bach measured in Munich (a small 
stream entering the Isar in Garching, north of Munich). Since we only consider the 
hydrological winter (November 1 to April 30), we have 2176 data points. 

First, we transform the margins to the unit interval applying the probability integral 
transform with the empirical marginal distribution functions. Then, we estimate the 
dependence structure using vine copula^: C\ 2 is estimated to be a Frank copula with a 
Kendall’s r o/ti 2 — 0.76, c 23 is a Frank copula with r 2 3 = 0.23 and ci 3;2 is a Gaussian 
copula with ri 3;2 = —0.18. Now we are interested in the resulting copula density of the 
maxima for one day (n = 2), one week (n = 14), one month (n = 60) and one winter 
(n = 362,). The respective contours (on the z-scale) are plotted in Figure 0 Similar 
to the examples from above we see that with increasing n the observed dependence 
structure tends to the independence copula (gray contours in the background). In case 
of the considered rivers this means that the maximal differences of the 12 hour average 
water levels over the entire winter are almost independent. 


3 Copula Density of Scaled Block-Maxima 


Examples 12.21 and 12.41 show that scaling of the block-maxima is necessary to achieve 
a limiting copula. These limiting copulas are called extreme value copulas and are 
ch aracte ri zed by m ax-sta bility. A recent introduction to extreme value copulas is given 
by Guden dorf and Segers ( 120101 ) . 

Since Hiisler and Reiss (j.1989) derived the scaling for the block-maxima of the mul¬ 
tivariate normal distribution with standard normally distributed margins X\,..., 
to a non-trivial extreme value copula, we use the same marginal scaling for the block- 
maxima MA”) on the z-scale given by 


Wj n) := b n [Mf - b n ) , 

where b n satisfies b n — n- p{b n ) for ip the standard normal density. Univariate extreme 

1 In order to assure that the necessary integrals are numerically tractable we had to exclude some 
pair-copula families (e.g. the t-copula). 
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for Wj G M. The marginal density of W ^' s is given by 


/^w = r rl K + ^rM^ + t 1 


for Wj G M, j = 1,..., d. Since is a strictly increasing transformation of M- n \ 

the copula of := ...., is the same as the one of M^. Therefore, 

using (11.31) we obtain the following expression for the joint distribution of 

F w(n) («; 1 ,..., w d ) = P (\v[ n) 

= C M( „) (V [b n + ^ (b n + ^ 
c ( 4+ ir)■ - ’*( 6 ” + ^ 


Similar arguments as in Corollary 11.21 can be used to express the joint density of 
in three dimensions for n > 3 as 


/ w w(Wl,W2,W3) = i^ n + P - ) 


+ ra(ra - 1)C (u 1 ,u 2 ,u 3 ) n 2 - SiC (ui, u 2 , u 3 ) d 23 C {u u u 2 , u 3 ) 


+ d 2 C (ui,u 2 , u 3 ) di 3 C (ui, u 2 , u 3 ) + asC (iti, u 2 , U 3 ) <9i 2 C (iti, u 2 , W 3 ) 

+ n(n — l)(?r-2)C , (-ui,rt2,ti3) n_3 aiC'(Mi,n 2 ,n3) (3-7) 

• a 2 C'(ui,'U2,ri 3 ) a 3 C'(Mi,u 2 ,'n 3 ) }, 


where u 3 T 


K + 7T ) for J = 1,2,3. 


K 


(n) 


According to Hiisler and Reiss ( 1989 1. besides scaling the maxima M-["\..., 
it is also necessary to change the correlation matrix E (n) = (pi,j(n)) 1<:i j <d of the 
underlying joint distribution of standard normal random variables X\ 1 ... .~X n , over 
whose i.i.d. copies X t ^,..., X ijd: i = 1,..., n, we take the maximum. The correlation 
matrices E (n) have to satisfy the following condition 


(1 - pij(n )) • log(n) -)• AC as n —>• 00 , (3.8) 

where \ t .j G (0,oo) are some constants for 1 < i,j < d, i ^ j and = 0 for 
i — 1,..., d. Since Pij(n) = Pjj(n), we also have = X 33 for 1 < i, j < d. Note that 
(13.81) implies that Pi,j{n ) —> 1 as n —> 00 . The limiting distribution H\ of the scaled 
maxima depends on A := (A i,j)i <i3<d . 


In the following we will examine the three-dimensional case and choose values for 
Ai 2 , A 13 , A 2 3 G (0, 00 ). For simplicity, we set 


Pi,j (n) 


1 - 


log (n) 


(3.9) 
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for 1 < i.j < 3, n e N, such that Equation 13.81 is always satisfied. However, for 
arbitrary it is not trivial to decide whether we obtain a valid correlation matrix 
through this particular choice of Pi.j(n) for any n G N. By construction the matrices 

/ 1 p 12 (n) pi 3 (n)\ 

E(n) = p 12 (n) 1 p23{n) 

\Pi3(n) p 23 (n) 1 / 

are symmetric and have ones on their diagonals. The only property we have to check 
is whether E (n) is positive definite. For this we only need to check if the determinant 
of each leading principal minor is positive. Since 1 > 0 and 1 — Pi 2 (^) 2 > 0 are trivially 
satisfied, the only real requirement is that 

det(E(n)) = 1 - pi 2 (n)pi 3 (n)p 23 (n) - pi 2 (n) 2 - p i3 (n) 2 - p 23 (n) 2 > 0. 


Using Equation 13.91 we obtain that det(E(n)) > 0 if and only if 


2 (^ 12^13 


+ A 2 2 A 


+ ^ 13 ^ 


) - (^12 + ^13 + ^23) > 


2A? 2 A 2 3 A 2 3 

log(n) 


(3.10) 


We denote the left-hand side of Equation 13.101 bv h( Af 2 , A 2 3 , A 23 ). Since the right-hand 
side of Equation 13.101 is always positive, it can only be satisfied if h(A 2 2 , A 2 3 , A 23 ) > 0. 
If h(A 2 2 , A 2 3 , A 23 ) > 0, then (13.101) is satisfied for all n E N with 


n > n 


exp 


2A? 2 A 2 3 A 


2 

23 



+ 1 


where |_ - J denotes the floor function. Table [I] shows the values of h and n* (if existing) 
for 10 different combinations of A12, A 33 , A23. 


# 

X 2 

A 12 

X 2 

A 13 

X 2 

A 23 

/i.(A 2 2 , A 2 3 , A 23 ) 

n* 

1 

1 

1 

1 

3 

2 

2 

2 

2 

2 

12 

4 

3 

1 

2 

3 

8 

5 

4 

0.5 

0.5 

0.5 

0.75 

2 

5 

0.3 

0.2 

0.1 

0.08 

2 

6 

0.2 

5 

0.75 

-15.8 

- 

7 

15 

20 

15 

800 

76880 

8 

100 

0.1 

20 

-6376.01 

- 

9 

1.05 

0.21 

0.84 

0.71 

2 

10 

4 

3 

3 

32 

10 


Table 1: Different A-combinations and the corresponding values of h and n*. 


Hii sler and Reiss! (1989) derived this scaling for multivariate normal distributions. 


Since we want to apply the scaling to vines, we need to transform the parameters of 
the normal distribution (correlations) to the parameters of the vine. 

Considering the vine structure from (12. pp we further assume that the pair-copulas 
are one-parametric. Having fixed Ai 2 ,Ai 3 ,A 23 such that h( A 2 2 ,A 2 3 ,A 23 ) > 0, we can 
perform the following procedure for n > n*: 


1. Calculate p i2 (n), p i3 (n) and p 2 3 ( 71 ) with the help of (13.91) . 
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2. Determine the corresponding partial correlation via 


Pl3;2(w) 


Pnjri) - Pi2{n)p 23 (n) 
a/1 - pi 2 {n) 2 y/l - p2z(n ) 2 


3. Translate the (partial) correlations pi 2 (n), P 23 (n) and Pi3 ;2 (t 7) into (partial) Kendall’s 
r values 712 ( 71 ), 723 ( 77 ) and Ti 3 ; 2 (t 7 ) using the relation for elliptical distributions 

r = — arcsin(p). 

7 r 


4. Determine the parameters # 12 ( 71 ), # 13 ( 77 ) and # 23 ( 77 ) of the pair copulas from the 
corresponding r values]^] 


Recall that pi 2 (7i) —* 1, piziji) —* 1 and p 2 3 ( 77 ) —* 1 as 77 —* 00 . Therefore, we also 
have 7 - 12 ( 77 ) — > 1, Ti 3 (?i) — > 1 and 723 ( 77 ) — > 1 as 77 —* 00 . However, the behavior of 
convergence of Pi3 ;2 (t 7) and hence Ti 3;2 (77) is not trivial. We use (13.911 to obtain 


Pl3;2 (ti) 



—y 


\2 

A 12 


+ A^ 3 - A 2 


13 


2AioA 


12 / '23 


as 77 —x 00 . Thus, 


7"i3;2 (ti) —X — arcsin 

7T 


12 


+ Ai 3 - A 2 


13 


2Ai 2 A 


12 / '23 


as 77 —x OO. 


For illustration, we will now take combinations 9 and 10 from Table [I] We show 
the (partial) correlations from Step 2 of the above procedure as well as the (partial) 
Kendall’s r values since they can be compared independently from the choice the 
respective pair-copulas. 



77 

P 12 ( 77 ) 

P 23 W 

Pl3;2 (Ti) 

7 - 12 ( 77 ) 

723 ( 77 ) 

7 - 1352 ( 77 ) 

Combination 9 

10 

0.54 

0.64 

0.87 

0.37 

0.44 

0.67 


50 

0.73 

0.79 

0.88 

0.52 

0.57 

0.69 


10 3 

0.85 

0.88 

0.89 

0.64 

0.68 

0.69 


OO 

1 

1 

0.89 

1 

1 

0.70 

Combination 10 

10 

-0.74 

-0.30 

-0.82 

-0.53 

-0.20 

-0.61 


50 

-0.02 

0.23 

0.25 

-0.01 

0.15 

0.16 


10 3 

0.42 

0.57 

0.44 

0.28 

0.38 

0.29 


OO 

1 

1 

0.58 

1 

1 

0.39 


Table 2: Overview over the (partial) correlations and (partial) Kendall’s r values for 
different n for combinations 9 (Af 2 = 1-05, Af 3 = 0.21, A 2 3 = 0.84) and 10 (Af 2 = 4, 
A 13 = 3, A 23 = 3). 

If we compare the values from Table |2] for combinations 9 and 10, it is eye-catching 
that the choice of A 12 , A 13 and A 23 has a crucial influence on the behavior of the (par¬ 
tial) correlations and the (partial) Kendall’s r values. In the first case the parameters 
are already relatively close to their limiting values for n = 10 3 , whereas in the second 

2 In the VineCopula package this transformation can be performed by the function BiCopTau2Par. 
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case they are still rather far from their limits for n = 10 3 . Further, we see that the 
limiting values of pi 3 - 2 (n) and ri 3 ; 2 (n) can be very different depending on the choice 
A 12 , A 13 and A 23 . 

Now we examine the behavior of the three-dimensional density of the scaled block- 
maxima / W (n> for increasing values of n. 


z 3 = <J> 1 (0.2) z 3 = <I> 1 (0.5) z 3 = <I> 1 (0.8) 





Figure 4: Two-dimensional slices of the density of the scaled block-maxima of a three- 
dimensional Clayton-vine (Af 2 = 1.05, Af 3 = 0.21, A 23 = 0.84). 


Example 3.1. First we look at a Clayton-vine and choose Xj 2 = 1-05, A 33 = 0.21 and 
A 23 = 0.84 (combination 9). The parameters and Kendall’s r values depend on the 
block-size n. 

Figure\4\ shows the density of the scaled block-maxima of the Clayton-vine for block- 
sizes n = 10, 50, 10 3 . Each row represents one block-size (and thus parameter set) and 
contains three contour plots corresponding to z^-values fixed to F~’ (n) (0.2), F~j n) (0.5) 

or F(^ n) ( 0 . 8 ), respectively. 
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z 3 = < J > 1 (0.2) z 3 = O '(0.5) z 3 = <jr'(0.8) 





Figure 5: Two-dimensional slices of the density of the scaled block-maxima of a three- 
dimensional Gaussian vine (A 32 = 4, A 33 = 3, A 23 = 3). 

Example 3.2. As a second example we choose a Gaussian vine with Xj 2 = 4, A 33 = 3 
and A| 3 = 3 (combination 10). 

Figured shows the density of the scaled block-maxima of the Gaussian vine for block- 
sizes n = 10, 50, 10 3 . Again, three contour plots corresponding to z^-values fixed to 
F~j n) (0.2), J n) (0.5) or F^} (n) (0.8), respectively, are displayed per row. The block-size 

and, Kendall’s r values are denoted on the left for each row. 


Conclusion 

In this chapter we showed that the copula density of the block-maxima of multivari¬ 
ate distributions can be expressed explicitly. For three-dimensional vine copulas we 
made use of the fact that we can compute their partial derivatives by one-dimensional 
integration, which makes the evaluation of the copula density for block-maxima numer¬ 
ically tractable. The advantage of our method is that we can use the entire sample for 
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estimation instead of reducing the sample size by taking the maximum over n observa¬ 
tions. Once we have estimated the underlying dependence structure we can derive the 
copula density of the block-maxima for any block-size (even larger than the original 
sample size). From the Clayton and the Gauss examples we have seen that without 
proper scaling the block-maxima do not approach a non-trivial limiting distribution 
for increasing block-size. 


Appendix 

Proof of Theorem 11.11 

In order to prove Theorem 11.11 we prove an auxiliary lemma from which Theorem 11.11 
follows as a corollary. 

Lemma 3.3. For k G { 1 ,..., d} and Uj G [0, 1 ] , j = 1 ..., d, we have 
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j=r [ 'P£S k j m&v 

Proof. We will prove this statement using induction. For k = 1 we have 
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The inductive step (k —> k + 1) proceeds as follows 


Qk+1 


du\ ■ ■ ■ du k+ i 
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du k+ i { du\ ■ ■ ■ du k 
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Applying the inductive assumption yields 
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We will consider the cases ”n > A:” and "n < /c” separately. We begin with Case 1 
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{n > k ): We have k An = k and hence 
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Now we use that fact that for k G {1,..., d — 1} and j G {1,... ,k An} we have 
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Applying Equation 13.111 yields 
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We perform an index shift in the first sum such that j + 1 is replaced by j and make 
use of the following two properties: 

(a) For all V G S^\ = {{{1,..., /}}} holds that {1} ^ V. 

(b) For all V G <S/j = {{{1} ,..., {Z}}} holds that {/} G V. 
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This results in 
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where we used the fact that k + 1 = (fc + 1) A n since n > k. This concludes the 
first case. Case 2 [n < k) is similar to the first one. The main difference is that 

C {v^ n ■,... ,«y n j = 1 for j = n, which was not possible before since j < k < n. 

Now k A n = n and therefore we obtain 
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where we applied Equation 13.111 in the second equality. In the third equality we per¬ 
formed an index shift in the Erst sum and used property (a) above. Since n < k we 
have n — (k + 1) A n. This concludes the second case and hence the proof of Lemma 

[3731 □ 

Having proved the auxiliary lemma we can now easily prove the statement from 
Theorem 11.11 


Proof of Theorem 11.11 Using Equation 11.31 we obtain 


c M («)( m i> •••,«<*) = 


3 d 


(Tl) • • • ) ^d) 


d d 




l/rix 
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du x • • • du d iV1 dm • • • du d L 

Now Theorem 11.11 follows directly from Lemma 13.31 for k = d. 

Proof of Theorem 12.11 

Proof. Expression 4) follows directly by the definition of c(ui,u 2 ,u 3 ) (see Equation 
[2761) . For expression 3(a) we can write 
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Expression 3(a) follows similarly. For expression 2(c) we have 
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Further, d 3 C(u\, 112 , 113 ) can be written in another way which we use for the calculation 
of the copula C(ui, u 2 , U3). In particular we have 
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Similarly one can derive expression 2(a). For the copula in expression 1 we have 
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Expression 2(b) follows by differentiating expression 1 with respect to v 2 . It remains 
to show expression 3(b). For this we differentiate 2(c) with respect to V\ . Therefore 
we have 
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Ci3;2(^112(l«l 1 ^ 2 ), C3|2(M 3 k 2 ))Cl 2 (ui, V 2 )c 23 (v 2 , U 3 )dv 2 . 
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